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(<-^ ' The Bethe-Salpeter equations for diquarks with the quantum numbers of the Nambu-Goldstone 

#^ , bosons are analyzed in the color-flavor locked phase of cold dense QCD with three quark flavors. 

The decay constants and the velocities of the Nambu-Goldstone bosons are calculated in the Pagels- 
Stokar approximation. It is also shown that, in contrast to the case of dense QCD with two flavors, 

Q, there are no massive radial excitations with quantum numbers of the Nambu-Goldstone bosons in 

the color-flavor locked phase. The role of the Meissner effect in the pairing dynamics of diquarks is 
^Nj ' explained. 
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C^ ■ I. INTRODUCTION 

f^ , Quark matter at high density is a color superconductor. It has been the subject of many studies for the last 

f^ • few years. This recent increased activity was initiated by the observation [|l|,y| that the color superconducting order 
^D ' parameter could be much larger than previously thought (for old studies, see Refs. [|[0]). Since then, many new 
(— I , results have appeared ||q-[20|]. 

^ ,' Because of asymptotic freedom, QCD becomes a weakly interacting theory at high densities. This allows one 

' to obtain some rigorous results for dense quark matter in the asymptotic limit. Of course, from the viewpoint of 

br^ phenomenology, it is most desirable to have a theory valid at intermediate densities that could be produced in heavy 
ri ion collisions or could exist in nature (for example, inside neutron or quark stars). This dilemma is partially resolved 

• • ' by studying predictions of the theory at high densities and, then, extending their validity as far as possible to the 
region of interest p-|ll||. Notice that all the heavy quark flavors could be safely omitted from the analysis when 
probing the quark matter at realistic intermediate densities. As a result, one arrives at a model of dense QCD with 
^ ■ either two ("up" and "down") or three ("up", "down" and "strange") flavors. 

In this paper, we deal with the problem of the diquarks with the quantum numbers of the Nambu-Goldstone 
(NO) bosons in cold dense QCD with three flavors. (A similar problem in the case of two flavors was considered in 
Refs. pl|,p3|.) The ground state of dense quark matter with three flavors is a color superconductor in the so called 
color-flavor locked (CFL) phase [O. It is remarkable that the chiral symmetry in such a phase is broken and most of 
quantum numbers of physical states coincide with those in the hadronic phase. It was tempting, therefore, to suggest 
that there might exist a continuity between the two phases |l^ . 

The dynamical symmetry breaking is caused by the celebrated Cooper instability in the pairing dynamics of quasi- 
particles around the Fermi surface. As a result, the original gauge symmetry SU(3)c and the global chiral symmetry 
SU{3)l X SU{3)r break down to the global diagonal SU{3)c+l+r subgroup |12|. Out of total sixteen (would be) 
NG bosons, eight are removed from the physical spectrum by the Higgs mechanism, providing masses to eight gluons. 
The other eight NG bosons show up as an octet [under the unbroken S'[/(3)c-i-l+_r] of physical particles. In addition, 
the global baryon number symmetry as well as the approximate U{1)a symmetry also get broken. As a result, an 
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extra NG boson and a pseudo-NG boson appear in the low energy spectrum. These particles are both singlets under 

SU{3)c+L+R- 

The low energy effective theory in the CFL phase was formulated in Refs. [p3|~p8[. Moreover, all parameters that 
define the low energy action of the (pseudo-) NG bosons were calculated in the limit of the asymptotically large 
chemical potential. The elegant method of Refs. pJ-Eq] is based on matching some vacuum properties (such as 
vacuum energy and gluon screening) of the effective and microscopic theories. While being rather efficient for some 
purposes, such a method is very limited when it comes to determining the spectrum of bound states other than NG 
bosons. The other approach used to study (diquark) bound states is based on the Bethe-Salpeter (BS) equation. In 
particular, the analysis of the BS equation seems to be the only way to test the conjecture of Ref. Q about the 
existence of an infinite tower of massive diquark states in the color superconducting phase of dense quark matter. In 



two flavor QCD, such a conjecture proved to be partially true ||2l],|2^. Here we generalize the method of Refs. ^1 22 
to the case of the CFL phase in QCD with three flavors 



This paper is organized as follows. In Sec. O, we describe the model and introduce the notation. Then, in Sec. Ill, 
we review the approach of the Schwinger-Dyson equation in the color superconducting phase of three flavor QCD. In 
Sec. |V|, we derive the Ward identities for the quark-gluon vertex functions, corresponding to the broken generators 
of global and gauge symmetries. We outline the general derivation of the Bethe-Salpeter equation and present its 
detailed analysis for the diquark NG bosons in Sec. ^ The decay constants and velocities of the NG bosons are 



calculated in Sec. VI. In Sec. VI], we discuss the problem of massive radial ex citat ions in the diquark channel with 



quantum numbers of the NG bosons. Finally, our conclusions are given in Sec. VIII 



II. MODEL AND NOTATION 

As we mentioned in the Introduction, the original SU{3)c x SU{3)l x SU{3)u symmetry of three flavor dense QCD 
breaks down to the global diagonal SU{3)c+l+r subgroup. The condensate in the CFL phase is given by the vacuum 
expectation value of the following diquark (antidiquark) field |12| : 

(0| {^n)';i' {^%)] |0) = ^,S!^S'^^ + n^S^Sl (1) 

where '^u and ^^ — C^^ are the Dirac spinor and its charge conjugate spinor, and C is a unitary matrix that 
satisfies C^^^^C — ~7j and C — — C^. The complex scalar functions ki and K2 are determined by dynamics. Here 
we explicitly display the flavor (i,j = 1,2,3) and color (a, 6 — 1,2,3) indices of the spinor fields. Notice, however, 
that the CFL phase mixes color and flavor representations. As a result, the notions of "color" and "flavor" become 
essentially indistinguishable in the vacuum. In passing, we also note that the order parameter in Eq. (0) is even under 
parity. 

By following Ref. |p4|, we introduce the color-flavor locked Weyl spinors [octets and singlets under SU{i)c+L and 
to replace the spinors of a deflnite color and flavor, 



S'[/(3)c+_R, respectively 



^A ^ i=7._(vl/g)/ (A^); , V3 = -^7'-(*g)/V, (3) 

~^A ^ -L^^^vl/g)/ (A-4); , ^ ^ _L^^(vl.g)/V, (5) 

where ^ = 1, . . . , 8, the sum over repeated indices is understood, and V± = (1 ± 7^)/2 are the left- and right-handed 
projectors. Tilde denotes the charge conjugate spinors. 

The order parameter in Eq. ffl) is recovered by assuming that the following (singlet under the locked symmetry) 
vacuum expectation values are non-zero: 

{0\i^m = -(O|0^|O) = -i (3^1 -f K^) , (6) 

(OIV^-^V^^IO) - -(O|<^^0^|O) = -\5^''>^2- (7) 



A specific vacuum alignment leads to a relation between ki and K2- For example, if the condensate were pure 
antitriplet in color and antitriplet in flavor (with respect to the original symmetries of the action) , it would imply that 
K2 = — Ki- Similarly, in the case of a sextet in color and sextet in flavor condensate, the relation would read K2 = ^i. 
It is known, however, that the true vacuum alignment is such that the antitriplct-antitriplet contribution dominates. 
At the same time, the sextet-sextet contribution is small but non-zero ||12 
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III. SCHWINGER-DYSON EQUATION 

In this section we briefly review the method of Schwinger-Dyson (SD) equation using our notation. This would also 
serve us as a convenient reference point when we discuss the Bethe-Salpeter (BS) equations in Sec. M. 
To start with, let us introduce the multi-component spinor, 



* 
^ 



A }, (8) 



built of the left-handed Majorana spinors, 



* = ^-hV', (9) 

^^^^^ + ^^^ (10) 

where A — 1, . . . , 8. Similarly, we could introduce the multi-spinors made of the right-handed Majorana fields, $ and 
$ . In our analysis, restricted only to the (hard dense loop improved) rainbow approximation, the left and right 
sectors of the theory completely decouple. Then, without loss of generality, it is sufficient to study the SD equation 
only in one of the sectors. 

The benefit of using the notation in Eq. (H) is that the inverse full propagator of quarks takes a very simple 
block-diagonal form, 



Si\p) = -I (^ + Ai7°7' + AiT'- +Air+), (12) 

S2\P) = -* (^ + ^7%' - A2V- - A2V+) . (13) 

Here Ai^2 = A^2^p ^ ^i2^p ^^^ ^1,2 = 7°Ai,27°' ^"^^ ^^^ "on-shell" projectors of quarks. 



,±_ Ml _u "'P \ ^_.,o. 



A- = -(^l±^), d = jH (14) 



are the same as in Ref. mM. We note that the gaps, Ai and A2, enter the propagators (12) and (03) with opposite 
signs. 

After expressing the standard bare vertex of QCD in terms of Majorana spinors (Q) and (|l^), we arrive at the 
following matrix form of the vertex: 

/ J_5^c^5 \ 

T^" = ^"^ (^ j^S^B^5 1 (rfAiJ^^S „ ^fABC) j , (15) 

where d'^^'^ and f^^'^ arc defined by the (anti-) commutation relations of the GcU-Mann matrices, 

[A^,A^] =2tf^^^X^. (17) 

With all the ingredients at hand, it is straightforward to derive the matrix form of the SD equation, 

G-\p) = G^\p) + 4^a. J ^^^^^G{q)T^%q,p)V,,{q~p). (18) 



Here ^^^ and T^^ are the bare and the full vertices, respectively. The gluon propagator in the SD equation is the 



same as in Ref. |15 , except that the Meissner effect should, in general, be taken into account. 



By inverting the expression in Eq. (|ll|), we obtain the following representation for the quark propagator: 



where 



. 7°bo + e+)-A+ . 7>o-4)-(^^)* 



+ ^ t'^'-f-^ KV. ^ ^y^-^-l^l KV^, (20) 

Po-(ep) -l^il Po-(ep) -l^il 

. 7°(P0 + e+) + A+ ^_^ . 7>o-4) + (A^)\ +^ 

+ , ^>°-S) + A.- ^^, . 7°bo + s) + (A,-)- 

^ pg-(ep-)^-|A,-p^^'^++%g-(e,-)2-|A,-p^-^- ('') 



The bare propagator in Eq. dlj) is similar but with zero values of the gaps A^ (n = 1,2). 

In the (hard dense loop improved) rainbow approximation, both vertices in the SD equation are bare. By making 
use of the propagator in Eq. dl^) and the vertex in Eq. (|l5|), we derive the set of two gap equations p5|,p^, 



Ar(p) = -f.«./(^rn^: 

A2(p) - ^Tra. j ^rV+ [252(g) - S^{q)] V+^V^^iq - p), 



2iq)v+rVf,4q-p), (22) 



(23) 



where we used the identities 



^ACD^BCD ^ 5^AB^ (24) 

jACDfBCD ^ 3^AB_ (25) 

An approximate analytical solution to this coupled system of gap equations was presented in Ref. |15| (a numerical 
solution was also given in Ref. IQ). Here we quote the final results | [l5| |, 

Ar(p4) ~ 2A^ipi) « 2A-3 3)(p4), (26) 

where 

A(3,3)(P4)«Ao, P4<Ao, (27a) 



A-g 3^(P4) « Ao sin I ^ ^ In A J , P4 > Aq, (27b) 

with A = 16(27r)3/2^/(3a^)^/2 ^j^^i 

16(2^)3/^ f 3^3/2 X 

As in the case of two flavor dense QCD, the issue of the overall constant factor in this expression is still unsettled. 
Some sources of possible corrections are discussed in Refs. |0,0|. In addition, we could argue (along the same lines as 
in Appendix B of Ref. [^) that there is also at least one non-perturbative correction that could modify the constant 
factor in the expression for the gap. 



IV. WARD IDENTITIES 

In order to preserve the gauge invariance in dense QCD, one has to make sure that Green functions satisfy the 
Ward identities. In this section, we consider the simplest Ward identities that relate the vertex functions and the 
quark propagators. In addition to establishing the longitudinal part of the full vertex function, these identities will 
play a very important role in our analysis of the BS equations for the NG bosons. 

In general, the structure of Ward identities in non-Abelian gauge theories (the Slavnov- Taylor identities) is much 
more complicated than in Abelian ones: they include contributions of the Faddeev-Popov ghosts. Fortunately in the 
(hard dense loop improved) ladder approximation, used in this paper, the situation simplifies. Indeed, since the direct 
interactions between gluons are neglected in this approximation, the Ward identities have an Abelian-like structure. 

To start with, let us rewrite the conserved currents (related to the baryon number and color symmetry) in terms 
of the Weyl spinors, defined in Eqs. (^ and (^0|). In this approximation, the left- and right-handed sectors of the 
theory decouple and the axial U{1)a charge is conserved. Therefore we can consider the (approximately) conserved 
currents in the two sectors separately. Here we give the details of the analysis for the left sector. The expressions for 
the singlet (related to the baryon number and the axial charge) and the octet (related to the color charge) channels 
read 

j^(.t) = ^oixh^.V+'^Dix) = h^{xh,j'^^{x) + hixh^^^'^ix), (29) 

Jp (a^) = \^D{xh^.X^V+^D{x) 

= i^^(x)7^ (d-^^^7' - */^''^) ^''(x) + ^ [^{xh,i'^^{x) + ^^{xh,j'^{x)] . (30) 

Similar expressions could be written in the right-handed sector too. Now, we are interested in the Ward identities that 
relate the quark-gluon vertices to the propagators of quarks. Therefore, let us introduce the following (non-amputated) 
vertex functions: 

T^{x,y) = (0|Tj40)*(x)*(2/)|0), (31a) 

r;t^(x,y) = {0\Tj,{Q)^^{x)^^{y)\0), (31b) 

r;^'^^(x, y) = {0\Tj^{Q)^^{x)^''{ym, (31c) 

Kfi^^y) = (o|rj-;f(o)*^(a:)§(2/)|o), (sid) 

rtf{x,y) = {0\Tj^{0)^ix)^^{y)\0). (31e) 

In order to derive the Ward identities, one needs to know the transformation properties of the quark fields. By making 
use of the transformation properties of the Dirac spinors, it is straightforward to derive the following infinitesimal 
baryon conservation symmetry transformations for the spinors of interest: 

5*-^ = icJ7^*^, (32) 

5*-^ = icJ*^7^ (33) 

(5* = iw7^*, (34) 

(5* = icJ*7^ (35) 

as well as the following color symmetry transformations: 

S^^ = -u,^ id^^'^r' -K z/-^^'^) vl/^ + ^7^^/, (36) 

2 V6 

<5^^ = ^c^^^^ (d^^^75 + z/^^^) + ^^7^, (37) 

2 V6 

,5* = ^7^*^, (38) 

V6 



%ui 



A 



b^ = ^^^7^. (39) 

V6 

Now, by making use of the current conservation as well as the definition of the vertices in Eq. (p^) , we straightforwardly 
derive the Ward identities for the non-amputated vertices: 



P'T^ {k + P,k) ^i [-f^Si (k) + Si{k + P) 7^] , 
PT;^^ (fc + P, fc) = i<5^^ [7^52 (fc) + ^2 (A: + P) 7'] , 
pppA.BC (^ + p^k) ^ '-d^''^ [7^52 (fc) + ^2 (fc + P) 7'] 

-^/^'"^[^2(fc)-52(fc + P)], 

P^r^f (fc + P, fc) = -^5^^ [7^5i (fc) + 52 (A: + P) 7^] , 
P^r^f (fc + P, fc) = -J=5^^ [7^52 (fc) + 5i (fc + P) 7^] , 



(40a) 
(40b) 



(40c) 
(40d) 
(40e) 



where 6*1 and S2 are the quark propagators. In the leading order approximation where the wave function renor- 
mahzation corrections are neglected, the explicit form of the momentum space propagators is given in Eqs. (EO) and 

(I). 

At this point, let us note that the use of the non-amputated vertices was crucial for the derivation of the Ward 
identities. Other than that, the non-amputated vertices are not very convenient to work with. In fact, it is the 
amputated rather than the non-amputated vertices that are usually used in the Feynman diagrams. Similarly, it is 
the amputated vertices that will appear in the BS equation in Sec. M. The formal definitions of the amputated vertices 
read 



r^ {k + p,k) = 5f 1 (fc + p) r^ (k + p, k) 5f 1 (fc) , 

r;^'' (fc + P,k) ^ S^' {k + P) r;t^ (fc + P, k) S^' (fc) , 

r;^'^^ {k + P,k)^ 5-1 (fc -t- P) Tf^^ (fc + P, k) S^' (k) , 

Tff (k + P,k)= s^' (k + P) r^;f (fc + p, k) s^' (k) , 

Kf (^ + P^k)= S^' [k + P) r^;f (fc + P, k) S^' (k) . 

It follows from Eq. (EO) that they satisfy the following identities of their own: 

P'T^ (fc + P,k) ^t [Sfi (fc + P) 7^ + 7^51-1 (fc)] , 
p^r;^^ (fc + P,k)= iS^^ [s^' {k + P) 7^ + 75^2-1 (fc)] , 

p,^ABC (fc + p^k) = '-d^^^ [S^' [k H- P) 7^ + 7'^2"' (fc)] 
-\f^''''[S^'{k + P)-S^Hk)\^ 
p^ri^;f (fc + P, A:) = -^5^^ [^2-1 (fc + P) 7^ + I'S^' (fc)] 

p^r^,;f(fc + p,fc) 



n/6 



6^''[S^'{k + P)^^+^^S^'{k)] 



(41a) 
(41b) 
(41c) 
(41d) 
(41e) 



(42a) 
(42b) 



(42c) 
(42d) 
(42e) 



In the rest of the paper, we are going to use these Ward identities a number of times. Because of a relatively 
simple structure of the inverse quark propagators, the last form of the identities for the amputated vertices will be 
particularly convenient. 

Notice that in the limit P ^ 0, we obtain 



p'^r^ (fc + p, fc)|p^„ = 2 (Ai(fc)p+ - Ai(fc)p_) , 



p^r;t''(fc + ^,fc)|p^o 



25 



AB 



A2(fc)P--A2(fc)P+ , 



p^r;^'^^(fc + p,fc)L „ =d^^^ 



IP^O 



p^r^'f(fc + p,fc) 



P^r^;f(A: + P,fc) 



1 






(A2(fc)P_-A2(fc)P+y, 

(Ai(fc) - A2(fc)) V+ + (A2(fc) - Ai(fc))P_ 
(Ai(fc) - A2(fc)) V+ + (A2(fc) - Ai(fc))P_ 



(43a) 
(43b) 
(43c) 
(43d) 

(43e) 



These expressions show that all five vertices have poles at P = 0. Such poles indicate the presence of NG bosons 
that correspond to broken continuous symmetries in the theory. 

Let us clarify the exact origin of the poles in all the ve rtice s introduced. It is straightforward to show that the 
pole contributions in the first two vertices, defined in Eqs. (31a) and ( 31b ), indicate the appearance of the NG boson 
related to breaking of the baryon number as well as the axial symmetry: in this approximation, the U{1)a charge 
is conserved and the axial symmetry is spontaneously broken. This boson is a sing let w it h res pect t o th e unbroken 
SU{3)l+r.+c subgroup. Similarly, one can see that the poles of the vertices in Eqs. (31c), (31d) and (31e) are due to 
an octet of the NG bosons related to breaking of the color symmetry (because of the locking, the chiral symmetry is 
also broken). Notice that considering the right-handed sector will double the number of the NG bosons. 
One could also establish that the explicit pole structure of the vertices should read 



T^ik + P,k)\p^^ 



■^AB 



■^A.BC 



ik + P,k)\ 



ik + P,k)\ 



r^;f(A + p,fc) 



r^;"" (fc + P, k) 



■ 2,ti 



where we introduced the notation P, 



(x) 



P^Q 

P^O 
2 , 



^^V(fc,o), 



-^, r,(fc,0), 

pMpin) 
P^^P- 

%f— d-^^^7ro(fc,0), 
p(^)p,^ 

pMpi^ 



pMpu 

pMpM 



nAB 



7ri(fc,0), 



6''"TV2ik,0), 



(44a) 
(44b) 
(44c) 
(44d) 
(44e) 



decay constants are defined by 

ipMpin)^ 
(0|j;f(0)|P;B)=*5^^FWi.W. 

The BS wave functions in the coordinate representation read 



{Pq, —c^P) (with Cx being the velocity of the appropriate NG boson). The 

(0|j,.(0)|P) 



and 



riix, y- P) = e-'^(^+«)/2^(x -y;P)^ (0|T*(x)#(y)|P), 
r,'(a:, y- P) ^ e-'^(-+'^)/2r,'(x - y; P) = \{Q\T^\x)^^{y)\P) , 



'!T^^^\x, y; P) = <5^^e-*-f'(^+^)/2^o(a; - y, P) ^ -d'^^°{0\T^^'{x)^^{y)\P; B), A, P = 1, 
{x, y; P) ^ <5^^e-^^(^+^)/27ri(x -y,P)^ {Q\T^''{x)^{y)\P; B), A, P = 1, . . . , 8, 



^1 ( 

7r^(^)(x, y; P) ^ 5^^ e-^^^^+y^'^ir^ix - y, P 



|r*(x)§'4(2;)|P;P), AP = 1,...,8, 



(45a) 
(45b) 

(46a) 
(46b) 

(47a) 

(47b) 
(47c) 



in the singlet and the octet channels, respectively. The momentum representation of the corresponding wave functions 
are given by the Fourier transforms of the translation invariant parts. 



V. BS EQUATIONS FOR NG BOSONS 



The bound states and resonances should reveal themselves through the appearance of poles in Green functions. To 
consider the problem of diquark bound states in cold dense QCD, one has to introduce four-point Green functions that 
describe the two particle scattering in the diquark channels of interest. The residues at the poles of the Green functions 
are related to the BS wave functions of the bound states. By starting from the (inhomogeneous) BS equations for 
the four-point Green functions, it is straightforward to derive the so-called homogeneous BS equations for the wave 
functions. 

The NG bosons in the problem at hand are the composite diquark states. Altogether, there are seventeen NG 
bosons and one pseudo-NG boson [the latter is related to breaking of the approximate U{1)a symmetry]. As was 



already mentioned in Sec. IV, there is no difference between the pseudo-NG boson and the real NG bosons in the 
ladder approximation. By combining the NG bosons in left and right sectors of the theory, we could construct the 
corresponding set of scalar and pseudo-scalar states. We note, however, that the scalar octet is unphysical because 
of the Higgs mechanism. The pseudoscalar octet is built of real NG bosons related to breaking of chiral symmetry. 
In addition, there is a scalar singlet and a pseudoscalar singlet. Both of them are physical, but the first is a real 
NG boson, while the other is a pseudo-NG boson. While considering the BS equation in the left-handed sector, we 
could reveal only half of the total eighteen NG bosons. As is clear, they should belong to a singlet and an octet 
representations of SU{3)l+r+c- 

We would like to note that though the scalars from the octet are removed from the physical spectrum by the Higgs 
mechanism, they exist in the theory as " ghosts" [g^ , and one cannot get rid of them completely, unless a unitary 
gauge is found. In fact these ghosts play an important role in getting rid of unphysical poles from on-shell scattering 
amplitudes Ig^l- Wc will use covariant gauges: because of the composite nature of the order parameter, it does not 
seem to be straightforward to define and to use the unitary gauge here. 

The nine Majorana quark fields in Eq. (||) form a singlet and an octet with respect to the color-flavor locked group 
of the vacuum. Clearly, we could construct all kinds of different diquarks out of these building blocks. If there is 
an attraction in the appropriate channel, two constituent singlets could give rise to a diquark singlet. Similarly, one 
singlet and one octet could produce an octet of diquarks. Finally, out of two quark octets, one could construct a 
whole new set of additional diquarks. Indeed, by making use of the representation theory, 

8 (g) 8 = 1 © 8^ © 8a © 10 © fO ® 27, (48) 

we might expect six more multiplets. As was shown in the previous section, however, only the singlet and the first 
of the octets (the symmetric one) are relevant for the description of the NG bosons. One could also guess that these 
singlet and octet are the most attractive channels in the color-flavor locked phase of dense QCD. For this reason, we 
restrict our study of the general bound state problem to the analysis of only these two channels [though one should 
keep in mind a possibility of a mixing between those two octets (see Sec. |VB| )]. 

In order to derive the BS equations, we use the method developed in Ref. |22|. To this end, we need to use the 
following effective Lagrangian: 

/2 — ^ iTf I J. I 0„,5 t A T) I A T") \ iT. I 1T..4 ^ „/ I 0„,5 AT) A T) I iT.^ 



-e// 



\^ [i) -K ^7%^ + AiP^ + Ai7'+) * + i*-'^ [i> + [i'fr' - A2-P- - A27'+) *' 



+ i*^(a;)A;^7^ (d^^C'75 - i/^^^) ^^ {x) + -^A^ (§(x)7^75^^(x) + ^^(x)7^75^(a;)) , (49) 

plus the right-handed contribution. This effective Lagrangian is a starting point in the derivation of the BS equations 
for the wave functions introduced in Eqs. (Eq) and (E^). While dealing with the multicomponcnt spinor defined in 
Eq. (g), it is rather natural to introduce the corresponding matrix form of the BS wave function. We use X['p\ P) as 
a generic notation for such a matrix wave function. 

In the (hard dense loop improved) ladder approximation, the matrix form of the BS equation reads 

G-^ [v + f ) ^(p; P)G-^ (^ ~ f ) ^ "^™^ / (0^^''^^'^' ^^^''''^t^^'i - V). (50) 



where T)r^^ {^^P) is the giuon propagator and 7"*^ is the bare quark-gluon vertex given in Eq. ( 15 ) . This approximation 
has the same status as the rainbow approximation in the SD equation. It assumes that the coupling constant is weak, 
and the leading perturbative expression for the kernel of the BS equation adequately represents the quark interactions. 

A. BS equation in the singlet channel 

Now, let us consider the BS equation in the singlet channel. As was explained above, it is sufficient to consider only 



the left-handed sector in this approximation. As we know from studying the pole structure of the vertices in Sec. IV, 
the singlet channel should contain at least one bound state, the diquark NG boson related to breaking of the baryon 
number [or, what is almost the same in our notation, to the approximate U{1)a symmetry]. 

We start the derivation of the BS equation from establishing the structure of the wave function. By incorporating 
the definition in Eq. (Eq), we obtain the following matrix form of the (non-amputated) BS wave function in the singlet 
channel: 



Now, we substitute this into the BS equation ( |50| ) and arrive at the following set of equations: 



In order to simplify the analysis, it is rather convenient to introduce the amputated BS wave functions, 

7y(p, P) = S^' (p + I] r7(p, P)S^' (^ " f ) ' (54a) 



V'ip, P) = S^' [p + I) rj'b, P)^^"! (^p - I) 



(54b) 



In addition, let us restrict ourselves to the analysis of the the BS wave functions in the limit P ^ 0. As one could 
check, such a limit is well defined and it is consistent with the on-shell condition for the NG bosons. The BS equations, 
then, read 



l^{p) = ^™, j ^j'^j^S2{qW{q)S2{qh'rV,,{q - p). (55) 

l^\p) = InasJ-^^^^ (j'S^iqMq)S^iqh' + ^j' S^iqW iq)S2iqh' " ^S^iqW S^iq)^ ^V^^q-p). (56) 



Now, we would like to get a solution to these equations. As in the case of two flavor QCD |21 2^, one could use the 
Ward identities to solve the problem. 



By comparing the Ward identities in Eqs. (43a) and ( 43b ) with the pole structure of the vertices in Eqs. (44a) and 



( [44q ) , we derive the required structure of the BS wave functions in the singlet channel 

vip) - ^ {s^Hph' + I's^^p)) = -^ (K,(j>)v+ - A^ip)r.) , (57) 

^'(P) = J(^ {S,\ph' + i'S-\p)) ^ -^ [A2{p)V- - A2ip)V+) . (58) 

It is straightforward to show that this is indeed the solution to the set of the BS equations ( |55| ) and (|56|), provided 
that Ai and A2 are the solutions to the gap equations (|2^) and (p3[). 



We use this solution in Sec. VI in order to derive the decay constants and velocities of the NG bosons. 



B. BS equation in the octet channel 

By following the approach of the previous subsection, let us also analyze the set of coupled BS equations in the 
octet channel. Again, it is clear from the pole structure of the vertices discussed in Sec. |^ that the this channel 
contains at least one solution which corresponds to the octet of diquark NG bosons. 

We start from establishing the general structure of the matrix BS wave function. By making use of the definition 
in Eq. (Wn), we obtain the following form of the (non-amputated) wave function in the channel of interest: 



X^(vP)-( ° S^^n2{p,P) \ 



(59) 



Notice, that here we added the "antisymmetric" octet cr(p, P) [see Eq. (H^)] to the general structure because this 
latter might in general mix with the NG octet. Of course, in the weakly coupled limit when the wave function 
renormalization is close to 1, the admixture of c r-octet is expected to be negligible. The reason is that, according 



to the (approximate) Ward identity in Eq. (43c), the NG boson does not contain the antisymmetric contribution 
proportional to f^^'~^- This is somewhat similar to the situation with the cr-singlet in two flavor QCD |23]. 

After substituting the wave function { p9{) into the BS equation (pO|), we arrive at the following set of equations: 



xYT^^Aq-p), (60) 

ST' (p + f ) ^2(P, P)S^' {P-\I ^ r"^ / W)^^' U^ii^^ Ph' + ;|^'^"('^' ^)^' - 3y|7'V(g, P) j 

x-rv^Aq-p), (61) 



^2'' (p + f ) ^o(p, P)s^' (p - f ) ^ ^ / (04^'' (2 Vl^' [^i('^' P) + ^2(9, i')] 7'^ - I'Mq, Ph' 

-37ro((Z, P) + 3 [a{q, P)r' - 7'V(<z, P)] j 7"2?m'>(9 - P)^ (62) 

^2"' (p + f ) '^(p, ^)52-' {p-\l =^l ^^' r^ [^1^*^' ^)^' " ^'^2^'^' ^)] + ^^°('^' ^)^' 

-57''7ro(g, P) + 5j'a(q, P)j' - 3a{q, P) j Y'D^.M - P)- (63) 
In deriving these equations, we used the following identities: 

^AA'B'j^BB'C'^CC'A' ^ _\^ABC ^ (g4) 

^AA'B'^BB'C'jCC'A' ^ J_jABC^ (gg) 

^AA'B'jBB'C'jCC'A' ^ _l^ABC ^ (gg) 

rAA'B' rBB'C' rCC'A' _ 3 rABC /gy\ 

Following the same approach as in the case of the singlet channel, we introduce the amputated wave functions by 

7ri(p, P) = S^' {p + 1^ 7ri(p, P)ST^ (^ - f ) ' (6^^) 

^2(P, P) = ^r' (^P + I) 7r2(p, P)52-i (^ " f ) ' (^^^) 

^o(p, P) - ^2"' [p + f ) 7ro(p, P)^^"^ (^ " f ) ' (^^^) 

a(p, P) - S^' [p + 1^ <t(p, P)52-i (P - f ) ' (^^^) 

and consider their BS equations in the limit of a vanishing total momentum of the diquark NG bosons, P ^ 0. 
Therefore, we arrive at the following equations: 



Mp) = l^^^j-^A^'^ \l'S^{q)TT2{q) + ^^'S2{q)M<l) + i\l\s2{q)a{q)\ S^iqh' l"" V ^^q - p) , 



(69) 



T2(P) = ^^"^/^^''^'^2(<Z) Ui(g)5i(q)7'' + -|7ro(<z)52(g)7''- 3^1^(9)^2(9)] 7'^2?p.(9 - P), (70) 

'ro(p) = ^ / (07^ [2^17^ [^2(g)7ri(<z)5i(g) + 5i(g)7r2(g)52(g)] 7' - 7^52(9)^0(9)^2(9)7' 

-352(9)7ro(9)52(9) +3 [52(9)a(rj)52(<z)7' -7'^2(9)<t(9)52(9)] \YV^,{q-p), (71) 
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~5^'S2{q)Mq)S2{q) + 5^'S2{q)a{q)S2iqh'-3S2{q)a{q)S2{q)]YT^,AQ-p)- 



(72) 



The appr oximate solution to this set of equations is obt aine d by c ompar ing t he Ward identities in Eqs. (43c), ( |43c]| ) 
and ( [43c| ) with the pole structure of the vertices in Eqs. ( [44c| ) , ( [44d| ) and (44e). In this way, we arrive at the following 
ansatz: 



7ro(p) 
7ri(p) = 



2F( 



^ {S^\ph' + 7'S-\p)) = -^ {A2{p)V^ - A2{p)V+) , 



V6F( 



{S2\ph'+l'S^\p)) = 



1 



7r2(p) 



' {S^\ph' + l'S2\p)) ' 



y6F('^) 



V6FM 



(Ai(p) - A2{p)) V+ + {A2{p) - Ai{p))V- 
(Ai(p) - A2(p)) V+ + (A2(p) - Ai(p))P_ 



a(p) = 0. 



(73) 
(74) 

(75) 
(76) 



It is straightforward to show that this is a solution to the set of the BS equations ( |69| ) - (jT^)- Notice, however, that 
the presented solution is approximate to the same extent as the quark propagators are. Indeed, if one takes into 
account the corrections due to the wave function renormalization of quarks, the expressions in Eqs. ([74|), (|75| ) and 
([Tq), determined by the Ward identities, would be modified. It is remarkable that no similar modifications would 
appear for the singlet NG boson considered in the previous subsection. This whole situation resembles quite a lot the 
analysis in two flavor dense QCD |23], where the NG doublets were free of any admixtures, while the NG singlet had 
a contribution from another singlet. All the arguments of Ref. |2|] in support of the selfconsistency of the leading 
order approximation apply without changes to the analysis here. 



VI. DECAY CONSTANTS OF NG BOSONS 



In this section, we derive the values of the decay constants for the NG b osons in t he Pagels-Stokar approximation 
pO| (for a review see Ref. Q). We start with the definitions in Eqs. (45a) and ( |45bD . It is straightforward, then, to 
derive the following exact expressions: 



/^ 9 



pM pM ^ __ 
^ 12 



12 



d^q 

(27r)4 

_dS_ 

(2^ 

d^q 

(2^)4 

d^q 



ti-{7M7'[8r7'(g,P)+r,(g,P)]} 

tr {7^7' [852 (q + P/2)v'{q, P)S2{q ~ P/2) + 5i(g + P/2)rKg, P)Sx{q - P/2)]} , (77) 

(57ro(<Z, P) + V67ri(q, P) + ^67^2(9, P)) 

^ ^hS2{q + P/2)7ro((?,P)52(9 - P/2) + \/652(g + P/2)7ri(g,P)5i(q - P/2) 

(78) 



tr 



tr 



7m7 



7m7 



(2^)4 

+V65i(q + P/2)^2(<Z, P)S2{q + P/2) 

In order to calculate the integrals on right hand side, one needs to know the explicit form of the BS wave functions at 
non-zero values of the total momentum P. In general, however, it is hard to derive them. Therefore, it is natural to 
consider the Pagels-Stokar approximation [pO| , pl[. Such an approximation uses the amputated wave functions at zero 
total momentum which are given in Eqs. (^7]), (|58|) and in Eqs. (73), (|7^) and (75) for t he singlet and the octet states, 



respectively. A simple calculation leads to the following result [see Eqs. (Al) and (A4) in Appendix A for details]: 



(pW)' 



^0 
clP 



27r2 



Po 

ip 
3 



7- 



lAri' 



|A2"I' 



247r2 
(21-81n2)//2 

72^1^ 



|Ari^ 



|A2-|2 



A2-|%jAri^ 

|A2-|2 



Po 

ip 
3 



Po 

ip 

3^ 



(79) 



(80) 
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where we used the relation between the gaps A]~ « 2A^ as in Eq. (^6|). 

Recall that here we consider the NG bosons from the left-handed sector. The expressions for the decay constants 
of the NG bosons from the right-handed sector are of course the same. In order to obtain the decay constants for the 
scalar and pseudoscalar NG bosons, one has to multiply the expressions for (F^''^)^ and (F'^'^^)^ in Eqs. (|79| ) and (jsy) 
by a factor of 2. 

We see that the decay constants of all NG bosons are of order /i, and their velocities are equal to 1/a/3. After 
taking into account the difference in definitions of the decay constants, we find that expressions d79) and (pO) agree 
with those derived in Ref. ||24[| P] (as well as with those in Ref. |g^, and up to a factor of 2 in (f")^ with those in 
Ref. PSJ), where a different approach was used. 

A few words should be said about Ref. |2^ , also using the method of the BS equation for studying the NG bosons in 
the CFL phase of cold dense QCD. The result for the decay constants obtained there is numerically different from our 
result as well as from the results of Refs. [ p4|p^j2^ . It seems the reason of that is because, while being qualitatively 
correct, the analysis of Ref. pSl misses some relevant quantitative details (for example, no distiction between the two 
types of quarks with different gaps seems to be made) . 

VII. ABSENCE OF MASSIVE RADIAL EXCITATIONS 

The important property of the quark pairing dynamics in cold dense QCD is the long range interaction mediated 
by the gluons of the magnetic type MM- Of course, the Meissner effect would eventually provide screening for the 
gluons in the far infrared region. It is known that such screening is negligible (at least in the leading order) for the 
dynamics of the gap formation R-nit]. The paring dynamics of the massive diquarks, however, is quite sensitive to 
the Meissner effect [|l]j22|. 

Let us now consider the difference between two flavor dense QCD and three flavor dense QCD. The former has the 
residual SU{2)c gauge symmetry in the vacuum, so that three out of total eight gluons do not feel the Meissner effect. 
In contrast, the gauge symmetry of the latter is completely broken (through the Higgs mechanism). This means, in 
particular, that all eight gluons in the CFL phase of three flavor QCD are affected by the Meissner screening. 

The natural question is whether the conjecture of Ref. jlj] about the existence of an infinite tower of massive 
excitations in the diquark channels with quantum numbers of the NG bosons could be generalized to three flavor 
QCD. We recall that this conjecture was derived from studying some unusual properties of the effective potential. 
Taken literally, the conjecture would imply the existence of an infinite tower of massive excitations for each of the 
five (would be) NG bosons in the case of QCD with two flavors. However, our detailed analysis ^,^ shows that, 
in fact, an infinite tower of the excitations occurs only in the color singlet channel. No radial excitations of the four 
other (would be) NG bosons (from a color doublet and antidoublet) appear. In order to reach this conclusion, one 
needs to take into account the effect of the Meissner screening for the pairing dynamics of diquark bound states. 

By following the same arguments as in the case of two flavor QCD [£l],g2|, one should distinguish between two 
classes of bound states, for which the role of the Meissner effect is very different. The first class consists of light 
bound states with the masses M <C |A|^|. The binding energy of these states is large (tightly bound states), and the 
Meissner effect is essentially irrelevant for their pairing dynamics. This point could be illustrated by the BS equations 
for the lightest diquarks, the massless NG bosons: in that case the most important region of momenta in the equations 
is given by \Aq\ < |fco| ^ \k\ ^ /x where the Meissner effect is negligible (see Sec. VII and Appendix C in Ref. [p2[). 

The second class includes quasiclassical states with the masses close to their threshold. Since the binding energy 
of the quasiclassical states is small, the quasiclassical part of the spectrum is almost completely determined by the 
behavior of the potential at large distances. In the particular case of cold dense QCD, the interaction between quarks 
is long ranged in the (imaginary) time direction and essentially short ranged in the spatial ones [pUl^fl- Because of that, 
the far infrared region, with |fco| < |Aq | ^ |fc|, is particularly important for the pairing dynamics of the quasiclassical 
diquark states and the Meissner effect is rather strong in that region. This implies that the inclusion of the Meissner 
effect is crucial for extracting the properties of the states from this second class (for more details see Sec. VII and 



Appendix C in Ref. l22|) 



Now, by repeating the arguments of Ref. [g2[, we conclude that, because of the Meissner effect, an infinite tower of 
(quasiclassical) massive diquarks does not appear in the CFL phase of three flavor QCD. This is directly related to 
the fact that all eight gluons in the model at hand are affected by the Meissner screening. This qualitative argument 



'^Our left-handed and right-handed NG fields correspond respectively to X and Y fields in the nonlinear realization of the 
SU{3)c X SU(3)l X SU{3)r symmetry of Ref. f^. 
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by itself does not prevent the possibility of a finite number of radial excitations in the spectrum. However, the same 
type of analysis as in Ref. |^ shows that no radial excitations of NG bosons appear at all. After the Meissner effect 
is taken into account, the interaction provided by gluons appears to be too weak to form even the lowest massive 
radial excitation of a NG boson. 

VIII. CONCLUSIONS 

In this paper we have studied the properties of diquark states with the quantum numbers of the NG bosons in the 
CFL phase of cold dense QCD with three quark flavors. We have derived the general Bethe-Salpeter equations in the 
singlet and the octet channels that include all the NG bosons. 

Our analytical analysis of the Bethe-Salpeter equations in the CFL phase shows that the theory contains one 
pseudoscalar octet of NG bosons, one singlet NG boson and one singlet pseudo-NG boson. This agrees with the 
previous results obtained in the effective theory approach of Refs. |p3|-p8| . We calculate the decay constants of the 
(pseudo-) NG bosons by using the Pagels-Stokar approximation, and the results agree with those of Ref. pj] , obtained 
by a different method. 

We also show that there are no spin zero massive states (with quantum numbers of the NG bosons) in the CFL 
phase of cold dense QCD. This conclusion is reached by studying a specific realization of the Meissner effect in three 
flavor QCD. In contrast to the case of dense QCD with two flavors, all eight gluons in the CFL phase become massive 
due to the Meissner effect. As a result, the one-gluon interaction becomes rather weak in the far infrared region which 
is responsible for the pairing dynamics of the (would be quasiclassical) massive radial excitations of the NG bosons 
p2| . We note, though, that our current conclusions state nothing about the bound diquark states of higher spins. 
The problem of the higher spin channels might be also interesting, but it is beyond the scope of this paper. 
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APPENDIX A: CALCULATION OF THE NG BOSON DECAY CONSTANTS 

In thi s ap pendix we exhibit in detail our calculation of the NG decay constants. We start from the definitions 
in Eqs. (|77| ) and (|78|). In the Pagels-Stokar approximation, the BS wave functions of the NG bosons are taken at 
vanishing total momentum P. The explicite form of the wave functions in the singlet channel is given in Eqs. (|57 
and (|5q). By making use of these expressions, we rewrite the right hand side of Eq. (^^ as follows: 

P('')i^(") ^-±-.J^tr{8j,S2{q)P[S,{q)+Y^S2iqh'] + J,S,{q) f' [S^iq) + ^'S,{qh']} + O {P') 

^J^pu [dq^( q,q,\ ( 8|A,-p \A^\' \ , 

-Fin)"" J (2.)^ 1^^^°^-- l#A[,| + (e,-)^ + |A-|f +[g| + (e,-)^ + |ArP]V ^^ 

-5^(5.oPo + Jp.)+0(P^), (Al) 

where we used the expansion of the quark propagators in powers of momentum P, 

Si^2{q + P/2) ~ Si,2{q) + 't;SiAq) PSi,2{q) + O (P^) , (A2) 



2 



as well as the following results for the Dirac traces: 



tr {-ff.A^j^A^V^) = 3^03^.0 - -prf • (A3) 
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In the case of the octet state, the calculation is similar. The explicite form of the corresponding BS wave functions is 



given in Eqs. (73), (|7J) and ([75[). By substituting them into Eq. ([78[), we arrive at 



[1 

[2; 
[s; 

[4] 
[5] 
[6] 
[7] 
[8] 
[9] 

[lo; 
[11 

[12 

[13 
[14 
[15 
[16 
[17 

[is; 

[19 
[20 

[21 
[22' 
[23^ 

[24' 
[25^ 

[26' 
[27 
[28' 
[29 
[30 
[31 



+ 27^ [S^{q) + 7^52 (g)7'] PS^{q) + 27^ [S^iq] + l'S,{q)j''] PS^ \ + O (P^) 



1 _ [ dq.d^q ( q,qA ( 5|A,-|^ , ^ 



9 1^09 1^0 



2 



+ 



|A~|^ A~(A^)* + A^(A^)* 

[ql + (e,-)2 + jA^-j^]' " [qi + {^^? + jArJ^] [g| + (e,-)^ + jA^'l^] 

M^ ^ Ar(A2-)- + A,'(Ar)- jAri^w i^\ , ,. 




M. Alford, K. Rajagopal, and F. Wilczek, Phys. Lett. B422, 247 (1998). 

R. Rapp, T. Schafer, E.V. Shuryak, and M. Velkovsky, Phys. Rev. Lett. 81, 53 (1998). 

B.C. Barrois, Nucl. Phys. B129, 390 (1977); S.C. Frautschi, in "Hadronic matter at extreme energy density", edited by 

N. Cabibbo and L. Sertorio (Plenum Press, 1980). 

D. Bailin and A. Love, Nucl. Phys. B190, 175 (1981); Nucl. Phys. B205, 119 (1982); Phys. Rep. 107, 325 (1984). 

R.D. Pisarski and D.H. Rischke, Phys. Rev. Lett. 83, 37 (1999). 

D.T. Son, Phys. Rev. D 59, 094019 (1999). 

D.K. Hong, V.A. Miransky, LA. Shovkovy, and L.C.R. Wijewardhana, Phys. Rev. D 61, 056001 (2000). 

T. Schafer and F. Wilczek, Phys. Rev. D 60, 114033 (1999). 

R.D. Pisarski and D.H. Rischke, Phys. Rev. D 61, 051501 (2000). 

S.D.H. Hsu and M. Schwetz, Nucl. Phys. B572, 211 (2000). 

W.E. Brown, J.T. Liu, and H.-C. Ren, Phys. Rev. D 61, 114012 (2000); ibid. 62, 054016 (2000). 

M. Alford, K. Rajagopal, and F. Wilczek, Nucl. Phys. B537, 443 (1999); M. Alford, J. Berges, and K. Rajagopal, Nucl. 

Phys. B558, 219 (1999). 

T. Schafer and F. Wilczek, Phys. Rev. Lett. 82, 3956 (1999); Phys. Rev. D 60, 074014 (1999). 

V.A. Mhansky, LA. Shovkovy, and L.C.R. Wijewardhana, Phys. Lett. B 468, 270 (1999). 

LA. Shovkovy and L.C.R. Wijewardhana, Phys. Lett. B 470, 189 (1999). 

T. Schafer, Nucl. Phys. B575, 269 (2000). 



S.R. Beane, P.F. Bedaque, and M.J. Savage, |nucl-th/00040i; . 

W.E. Brown, J.T. Liu, and H.-C. Ren, Phys. Rev. D 62, 054013 (2000). 

D.H. Rischke, Phys. Rev. D 62, 034007 (2000); ihid. 62, 054017 (2000). 

G.W. Carter and D. Diakonov, Nucl. Phys. B582, 571 (2000); R. Casalbuoni, Z. Duan, and F. Sannino, Phys. Rev. D 62, 

094004 (2000). 



V.A. Miransky, LA. Shovkovy, and L.C.R. Wijewardhana, liep-ph/0003327 



V.A. Mhansky, LA. Shovkovy, and L.C.R. Wijewardhana, Phys. Rev. D 62, 085025 (2000). 



R. Casalbuoni and R. Gatto, Phys. Lett. B 464, 111 (1999); bep-ph/9911223 ; D.K. Hong, M. Rho, and I. Zahed, Phys. 

Lett. B 468, 261 (1999). 

D.T. Son and M.A. Stephanov, Phys. Rev. D 61, 074012 (2000); ibid. 62, 059902(E) (2000). 

M. Rho, A. Whzba, and I. Zahed, Phys. Lett. B 473, 126 (2000); M. Rho, E. Shuryak, A. Wirzba, and I. Zahed, Nucl. 

Phys. A676, 273 (2000). 

D.K. Hong, T. Lee, and D.-P. Min, Phys. Lett. B 477, 137 (2000); C. Manuel and M.G.H. Tytgat, ibid. 479, 190 (2000). 

K. Zarembo, Phys. Rev. D 62, 054003 (2000). 

S.R. Beane, P.F. Bedaque, and M.J. Savage, Phys. Lett. B 483, 131 (2000). 

R. Jackiw and K. Johnson, Phys. Rev. D 8, 2386 (1973); J.M. CornwaU and R.E. Norton, ibid. 8, 3338 (1973). 

H. Pagels and S. Stokar, Phys. Rev. D 20, 2947 (1979). 

V.A. Miransky, Dynamical Symmetry Breaking in Quantum Field Theories (World Scientific, Singapore, 1993). 

14 



